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Abstract 



C\ • This paper presents recent results concerning the existence and qualitative properties of 

D ' travelling wave solutions to the Gross-Pitaevskii equation posed on the whole space R*. 

Unlike the defocusing nonlinear Schrodinger equations with null condition at infinity, the 
fSJ | presence of non-zero conditions at infinity yields a rather rich and delicate dynamics. We 

C^l i focus on the case N — 2 and N = 3, and also briefly review some classical results on the 

one-dimensional case. The works we survey provide rigorous justifications to the impressive 
p I i series of results which Jones, Putterman and Roberts [46, 45] established by formal and 

numerical arguments. 

I Introduction 

The Gross-Pitaevskii equation 

id t ^ = A^ + ^(l- |^| 2 ) onlR^ xR, (GP) 

""xl" ! 

appears as a relevant model in various areas of physics: nonlinear optics, fluid mechanics, Bose- 
Einstein condensation... (see, for instance, [57, 40, 44, 46, 45, 48, 12]). It corresponds to a version 
of the defocusing nonlinear Schrodinger equations. If one considers finite energy solutions, as 
■ suggested by the formal conservation of the energy (see (1.1) below), then \& should not vanish 

OS . at infinity, but instead \^(x, •)[ should in some sense tend to 1 when \x\ — > +oo. As a matter of 

fact, this condition ensures that (GP) has a non-trivial dynamics, contrary to the case of null 
condition at infinity, where the dynamics is expected to be trivial (dispersion, scattering...). 

In nonlinear optics, (GP) classically appears in the context of optical dark solitons, that is 
localized nonlinear waves (or "holes") which exist on a stable continuous wave background (see 
the review paper by Kivshar and Luther-Davies [48] or paper [49], where different scenarios for 
the transverse instability of one-dimensional black solitons are proposed in the context of (GP)). 
The boundary condition \^>(x, •)[ — ► 1 is due to this non-zero background. 

The Gross-Pitaevskii equation has also been intensively used as a model for superfluid Helium 

II and for Bose-Einstein condensation. We refer to [40, 21], and to the survey by Berloff [7]. 

At least on a formal level, the Gross-Pitaevskii equation is hamiltonian. The conserved 
Hamiltonian is a Ginzburg-Landau energy, namely 

E^) = \! \V*\ 2 + ][ (1-|*| 2 ) 2 ^/ eW . (1.1) 
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Similarly, the momentum 

P(*) = \ f (iV* (1.2) 

is formally conserved. We will denote by p, the first component of P, which is hence a scalar. 
Another quantity which is formally conserved by the flow is 



m («r) = - f (\^\ 2 - 1 



If does not vanish, one may write 

* = y^expzy?. 
This leads to the hydrodynamic form of the equation 

( d t p + div(pv) = 0, 

\ p(d t v + v.Vt/) + Vp 2 = pV - ^f) , (L3) 

where 

v = -2V<p. 

If one neglects the right-hand side of the second equation, which is often referred to as the 
quantum pressure, system (1.3) is similar to the Euler equation for a compressible fluid, with 
pressure law p(p) = p 2 . In particular, the speed of sound waves near the constant solution 
^ = 1, that is p = 1 and v = 0, is given by 

c s = y/2. 



When N = 1, (GP) is integrable by the inverse scattering method, and it has been formally 
analyzed (with the non-trivial boundary condition at infinity) in [61]. The (local and global) 
Cauchy problem associated to (GP) with non-zero condition at infinity has been recently thor- 
oughly investigated (see the survey by Gerard [30] in this volume). The possible existence 
of non-trivial travelling wave solutions highly suggests that the long time behaviour of global 
unsteady solutions should be non-trivial (as it would be in the presence of zero condition at 
infinity). 

Travelling waves are solutions to (GP) of the form 

ty(x, t) = u(xi - ct, x±), x± = (x 2 , • • • , xn)- 

Here, the parameter c£M corresponds to the speed of the travelling waves (we may restrict to 
the case c > using complex conjugation). The equation for the profile u is given by 

icdm + Au + u(l - \u\ 2 ) = 0. (TWc) 



Travelling waves of finite energy play an important role in the dynamics of the Gross- 
Pitaevskii equation. Whereas the one-dimensional case can be integrated explicitly, solutions 
are not explicitly known in higher dimensions. 

The first attempts to prove existence of travelling waves and to study their properties were 
based on formal expansions and numerics. They were mostly performed by physicists. The three- 
dimensional case was considered by Iordanskii and Smirnov [44], whereas Jones, Putterman and 
Roberts [46, 45] developed a thorough analysis in dimensions two and three: they found a branch 
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of solutions with speeds c covering the full subsonic range (0, \/2), and they conjectured the non- 
existence of travelling waves for supersonic speeds. Moreover, they showed the connections with 
solitary waves for the Kadomtsev-Petviashvili equation (KP I). 

More recently, the program of Jones, Putterman and Roberts [46, 45] has been provided with 
rigorous mathematical proofs. One of the purposes of these notes is to review some of these 
mathematical progresses and to stress a number of open problems. 

The next section provides a rather extensive survey of the one-dimensional case. We review 
classical results and also prove some new ones. Sections 3 and 4 are devoted to recent progresses 
on the rigorous mathematical proofs of the program of Jones, Putterman and Roberts [46, 45]: 
qualitative properties and existence of travelling waves in higher dimensions. Section 5 presents 
some related problems. In particular, we establish new results on the nonlinear Schrodinger flow 
around an obstacle, modelled by a potential. 



2 The one-dimensional case 
2.1 The integration of equation (TWc) 

In the one-dimensional case, equation (TWc) is entirely integrable using standard arguments 
from ordinary differential equation theory. Hence, it is possible to classify all the travelling 
waves with finite energy according to their speed c. 

Theorem 1 ([66, 54, 36, 29, 25, 7]). Assume N = 1 and c > 0, and let v be a solution of finite 
energy to (TWc). 

i) If c> v2, v is a constant of modulus one. 

ii) If0<c< y/2, up to a multiplication by a constant of modulus one and a translation, v is 
either identically equal to 1, or to 



2-C 2 , fy/2 



c 2 \ . c 



v(x) = t) c (x) = \l — - — th( x)-i-j=. (2.4) 



Moreover, D c may be written as 
c (x) = 



2-c 2 ./ /V2^ , /V2 



1 -. r- exp i arctan th I x I , (2-5) 



c 2 \ \ 7T 
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if < C < y/2. 

For sake of completeness, we briefly recall the proof of Theorem 1. 
Proof. Denoting v = v± + iv2, equation (TWc) becomes 

v'i-cv' 2 + v l (l-vf-vl) = 0, (2.6) 

t4' + cui+«2(l-u?-u£) = 0. (2.7) 

Multiplying (2.6) by v 2 and (2.7) by v\ yields 

M - v 2 v' 1 )' = y, (2.8) 
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where 77 = 1 — \v\ 2 . Since v is a finite energy solution to (TWc), it is smooth and bounded on R. 
Moreover, v'(x) and rj(x) tend to 0, as x — > ±00. In particular, the integration of (2.8) leads to 

v x v 2 - v 2 v[ = -77. (2.9) 

Multiplying (2.6) by v[, and (2.7) by v' 2 , and integrating the resulting equation, we also deduce 
that 

W\ 2 = y. (2-10) 

Computing 

rj" = -2\v'\ 2 - 2(v 1 v'{ + V2V2) = -2\v'\ 2 - 2c{v x v 2 - v 2 v' l ) + 2rj - 2-q 2 , 

where we invoke (TWc) for the second identity, we finally obtain from (2.9) and (2.10), 

77" + (c 2 - 2)77 + 3r] 2 = 0. (2.11) 

Multiplying (2.11) by 77' and integrating, it follows that 

rj 2 + (c 2 - 2)r} 2 + 2rf = 0. (2.12) 

As previously mentioned, 77 is a smooth, bounded function on ]R, such that r/(x) — > 0, as x — > ±00. 
In case it is not identically equal to on R, we can assume, up to a translation, that 

|t7(0)| =max{\r)(x)\,x el} > 0, (2.13) 

so that 

r/(0) = 0. (2.14) 

Using (2.12), we are led to 

n(o) = ' 2 -^-. (2.15) 

In particular, if c = \^2, 77(0) = 0, which gives a contradiction with (2.13). Hence, r\ is identically 
equal to if c = \[2. On the other hand, if c > \/2, then, by (2.12), 

(c 2 - 2 + 2t])t] 2 = -7]' 2 < 0. 

so that, by (2.15) and the continuity of 77, 

2 - c 2 

V(x) < < 0, Vx6 R. 

This gives a contradiction with the fact that 77 G L 2 (R). Therefore, 77 is identically equal to if 
c > \[2. Finally, if < c < \^2, it follows from (2.14), (2.15), and Cauchy-Lipschitz's theorem 
that there exists a unique local solution 77 to (2.11)-(2.14)-(2.15) on some neighbourhood of the 
point 0. One can check that this solution is given by 



rj{x) = 2 * s , 
2ch 2 (^pV 



so that it can be extended to be the unique, up to some translation, non-constant solution to 
(2.11) in L 2 (R). In conclusion, either 77 is identically equal to (in particular, when c > \/2), 
either it is given, up to a translation, by (2.16). 
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When c 7^ 0, it follows that \v\ = yj\ — rj does not vanish on R, so that we can construct a 
lifting ip of v, that is a smooth function (p such that 



v = \v\ expi(/? = y/l — r) exp 

Invoking (2.9), we are led to 



* 2-27?' v ; 



so that = (po is identically constant when rj = 0, whereas </? is given by 



V2^2 



c 



.2 



(^(x) = tpo + arctan I th^ x 

when rj is given by (2.16). This completes the proofs of assertion i), and of identity (2.5). 
Identity (2.4) follows from standard trigonometric identities, in particular, from the fact that 

cos ( arctan(x)) = , and sin ( arctan(x)) 



Vl + x 2 v ' y/l + x 2 

When c = 0, either r\ is identically equal to 0, so that, by (2.10), v is a constant of modulus 
one, or r\ is given, up to some translation, by 

r\(x) 



ch 



X 

\V2 



In this case, v(0) = 0, so that, by (2.10), we may assume, up to the multiplication by a constant 
of modulus one, that 

</(0) - 

Hence, by Cauchy-Lipschitz's theorem, there exists a unique local solution to (TWc), with c = 0, 
on some neighbourhood of the point 0. One can check that this solution is given by 

t)o(z) =th(-^=), 

so that, up to the invariances, 

v(x) = Oo(x) =th(-^=). 
This completes the proof of Theorem 1. □ 

Let us next draw some remarks concerning qualitative properties of non-constant finite energy 
travelling waves, which we deduce from formulae (2.4) and (2.5). 

The first observation is that the one-dimensional non-constant travelling waves form a smooth 
branch of subsonic solutions to (TWc). As seen in Theorem 1, there exist neither sonic, nor 
supersonic non-constant travelling waves. Moreover, formula (2.4) yields the spatial asymptotics 
of the non-constant solutions to (TWc). Notice in particular that 



* 0^°° = ±\l - i-^?=, as x — > ±00. 



Hence, t> c (x) converges to a constant of modulus one, as x — > ±00. This constant is different 
for the limit x — ► —00, and the limit x — > +00. In contrast, the solution we will find in the 
higher dimensional case will have a limit at infinity which is independent on the direction (see 
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Theorem 9). Notice also that the function d c — 0^°° has exponential decay, whereas the decay 
is algebraic in higher dimensions (see also Theorem 9). This fact is related to the properties of 
the kernel of the linear part of the equation. The Fourier transform of the kernel K associated 
to (TWc) (more precisely, to (2.11)) is given by 



2 - c 2 + £ 2 ' 

hence has some exponential decay, when < c < 

A second observation is, in view of identity (2.5), that D c (x) does not vanish, unless c = 0. In 
the case c = 0, the non-constant solution of finite energy to (TWc) is given, up to the invariances, 
by 

D (aO = th(-^=), 

which vanishes at x = 0. This stationary solution, which is real-valued, is often termed a kink 
solution. In contrast with the one-dimensional case, it can be shown in higher dimensions, using 
Pohozaev's identity, that the only finite energy solutions with c = are constants (see [11]). 
The energy of the kink solution is given by 

The kink solution has the following remarkable minimization property. 
Lemma 1. We have 

E{t> ) =inf \E(v), v e Hl c (R), influx)! = o). 
In particular, if E(v) < ^j^, then, 

inf \v (x)\ > 0. 

Proof. We consider a minimizing sequence (f n )neN for the minimization problem 

So = mfy^°°e(v),v G fl£ c ([0,+oo)),i;(0) = oj, 

which is well-defined by Sobolev's embedding theorem. We notice that the functions v' n are uni- 
formly bounded in L 2 ([0, +oo)), and that v n (0) = 0. Hence, by Rellich's compactness theorem, 
there exists some function u G i7j^ c ([0, +oo)), with u(0) = 0, such that, up to a subsequence, 

v' n — v v! in L 2 ([0, +oo)), and v n — > u in LJ^ c ([0, +oo)), as n — > +oo. 

By Fatou's lemma, we are led to 

so that the infimum £$ is achieved by the function u. In particular, the solution u is critical for 
the Ginzburg-Landau energy, i.e. solves 

u" + u(l- \u\ 2 ) = 0. 
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Integrating this equation yields u(x) = t>o(x) = th^-^=j (see the proof of Theorem 1), so that 

f+°° sj^ 
£o = y o e(t» ) = - g-. 

Next consider a map v £ fl^ (R), with finite Ginzburg-Landau energy, and which vanishes at 
some point xo- In view of the invariance by translation, we may assume that xq = 0, whereas 
the minimality of Do yields 

e(u) > — , 
o 6 

and the same inequality holds for the energy on (— oo,0], so that 

E(v)>^- = E(x>o). 

The proof of Lemma 1 follows. □ 
Remark 1. When E(v) < ^p, v has no zero, so that we may write 

v = £>exp iip. 

In particular, the energy density of v can be expressed as 

e(v) = \{{ e '? + e\v'?)+\^-8 2 )\ (2.18) 



2 V u ' " XT J J 4 
whereas, for the scalar momentum, we have 

(i v ', v ) = -£>y. (2.19) 
2.2 Energy and momentum of travelling wave solutions 

A general principle states that travelling wave solutions are critical points of the energy, keeping 
the momentum fixed. A natural question is therefore to determine whether this principle applies 
to the solutions t> c . However, the precise definition of the momentum raises a serious difficulty. 

In the context of nonlinear Schrodinger equations, the momentum of solutions should be 
defined as 



This quantity has a physical meaning and at least formally, describes the evolution center of 
mass in the (GP) dynamics. This fact is shown for instance by the relation 



d fl 

dt\2 



j x(\V(x,t)\ 2 -l)dx] =2P(*(i)). (2.20) 



The momentum P(t> c ) is well-defined in view of the exponential decay of the quantity (io' cl d c ), 
and is equal to 

P(*c) = — • 

Notice that, if c ^ 0, then we may write D c = g c expi(p c , so that we have 



P(f c ) = -\ [ fid- 
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In the variational context we have in mind, we need to extend this definition to a larger class of 
maps v. S We restrict ourselves to the case the case E(v) < ^p, since we wish to establish a 
minimality property for the solution d c , which is known to have an energy satisfying E(t) c ) < 2 4p. 

Assuming E(v) < ^p, it follows from Lemma 1 that v has no zero and the continuity of v 
combined with the finiteness of its energy yield 

5 = inf \v(x)\ > 0, 
so that we may write v = gexpiip. However, the quantity 

z JR z JR z JR 

is only well-defined if we are able to give a meaning to the last integral L <p', since the first one 
may be bounded by Cauchy-Schwarz's inequality, 



/ (i - ev 

JR 



The integral L ip' is well-defined in particular if the limits <^(±oo) = lim <p(x) do exist, in 

x— >±oo 

which case we have Ly' = ip(+oo) — (p(—oo). This leads us to consider the set 

Z (R) =\v€ HlJR), s.t. E(v) < and v ±00 = lim v(x) exist). 

I 3 x—>±oo ) 

In view of the above discussion, the momentum P(v ) is well-defined for any v G Zq(M.). However, 
we have 

Lemma 2. Let p > 0. We have 

inf Ie{v), v £ Z (R), s.t. P(v) = p } = 0. 

Proof. Let n > 1, and consider the map v n = expiip n G H^ oc (M) defined by 

{0, if x < — n, 
_ p x±n ) if _ n < a; < nj 
— 2p, if x > re. 

Since 

P 2 

E(v n ) = - >0, as n -> +oo, (2.21) 

n 

v n (x) —* 1, as x — oo, and v n (x) — > exp(— 2ip), as x — > +oo, the map v n belongs to Zq(M.), for 
any n sufficiently large. Moreover, we compute 

P(v n ) = p, 

so that Lemma 2 follows from (2.21) taking the limit n — ► +oo. □ 

We do not actually know if the space Zq(M), and the limits at infinity for maps in that space 
are preserved by the (GP) flow. However, this is the case if one imposes more regularity on the 
initial condition. Consider, for instance, the space 

Z 2 (R) = Zq(R), s.t. U|( ,±oo) G {v ±oo } + H 2 (0,±oo)Y 

We have 
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Lemma 3. Assume 6 ^(K), anc ^ ^ ^ ^ e ^ e solution to (GP) wrai/i initial datum 'J/q- 
Then, *&(t) belongs to Z<i QR) /or any i £ R, and Vl/±oo(i) = ^o,±oo- 

We will not provide a proof of Lemma 3 in this survey. 

Lemmas 2 and 3 lead us to introduce a notion of renormalized momentum p, whose physical 
meaning seems less obvious, but whose variational properties are more relevant for our study. 
We set, for any function v £ H^ oc (R) such that E[v) < ^p, 

P(v) = l I (1 - eV, (2.22) 



2 _ 

where we use the polar form v = gexpiip of v. Notice that, by construction, 

P(v) = p( v ) - - (</?(+oo) - <^(-oo)) , 

when the limits <^(±oo) exist. We next compute the values of the energy and renormalized 
momentum. 

Proposition 1. Assume N = 1 and < c < \[2, and let t> c be the non-constant solution of 
equation (TWc) given by (2.4). Then, the energy of tt c is equal to 

(2 - r 2 ^ 

E(X> C ) = { 3 ' , (2.23) 
whereas its renormalized momentum is given by 



p{* c ) = -- arctan \q=^) ~ ^T^- ( 2 - 24 ) 



Remark 2. It should also be noticed that 

,2 c 2 - 2 

2ch 2 (^x) 

and that this quantity is integrable, so that the mass of d c is equal to 



m(t> c ) = -yjl- c 2 , 

whereas its mass center is given by 

\ [ x(\* c (x)\ 2 -l)dx = 0. 
Remark 3. In view of identity (2.24), it is tempting to pass to the limit c — > 0, and state that 

p(do) = \ (2.25) 
As a matter of fact, this identity might be recovered using the singular lifting 

tJo(x) = th ("^7|) ex P (^(sig 11 ^) - l)|)> 

so that the derivative of the phase ipQ is equal to tt5q, an d formula (2.22), in the sense of measures, 
yields (2.25). 
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Proof. It follows from (2.4) and (2.10) that 
2|*' c (*)| 2 = (l-| Dc (x)| 2 ) 2 = 

so that 



c 2 \2/ / V2 - c 2 \2 
1 1-th x 



E(t> c ) 



(2 - c 2 )S 



l-th(u) 2 du 



(2-c 2 )§ 



Using definition (2.22) and identity (2.17), we next compute 

4 ii 1 



c 



»7c + 



where 7] c = 1 — |t) c | 2 . By (2.4), we are led to 

P(t>c) 



cV2 - c 2 



+oo 



:V2 - c 2 ^2" 



1 -th(it) 2 )<fot + 2 
^ 2 /" 1 dv 



+oo 



1 - W 



1 - th(u) 



c 2 + (2 - c 2 )th(<u) : 



;du 



1 + ^1 



dv 



V2" 



arctan 



V2 



which gives formula (2.24). 



□ 



Invoking formula (2.24), we notice that the function c i— » p(t) c ) is smooth, decreasing, and 
satisfies 

±(p(p c ))=-V2=#. (2.26) 

Hence, it performs a diffeomorphism from (0, y/2) on (0,5), so that we can express -£7(t> c ) = 
£ (p(t c )) as a function of p(t> c ) to obtain the following graph. 



E 




The curve E = £{p) is smooth, increasing and strictly concave, and lies below the line 
E = y/2p. Each point of the curve represents a non-constant solution D c to (TWc) of energy 
E(y> c ) and scalar momentum p(v c ). The speed of the solution D c (and as a result, its position 
on the curve) is given by the slope of the curve. Indeed, it follows from (2.23) and (2.26) that 

|(*))^K))(^M)"- 
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2.3 The variational interpretation of equation (TWc) 

It is a rather general principle that travelling waves could be obtained minimizing the Ginzburg- 
Landau energy keeping the momentum fixed. In view of Lemma 2, the appropriate notion of 
momentum is given by formula (2.22). In this context, we consider the set 

X p = [v £ fl£ c (R),s.t. E(v) < Hp, and p(v) = p}, 

for any p > 0, and the minimization problem 

E ain (p) = -wf{E(v),v€X p }. (V p ) 

For any < p < ^, let c = c(p) be the only speed c such that 

p = arctan ( ■ =] -J 2 — c 2 . 

2 V^2^/ 2 V 

We have 

Theorem 2. Ze£ < p < \. Then, Xp is not empty, and we have 

E(t) cip) ) =inf {e(v),v £ X p }, 

that is (Vp) is achieved, and the only minimizer, up to invariances, is t) c ( p ). 

Remark 4. For p = 0, the solutions to the minimization problem (Vp) are the constants of 
modulus one, for which 

E min (0) = 0. 

Remark 5. Notice that the set Xp is empty when p > f . It is formally tempting to assert that 

#min(P) — — — , 

for any p > -|. Indeed, any map which vanishes at some point has an energy larger than 3^, 
and we can construct minimizing sequences with the appropriate limiting energy. For instance, 
we may consider the functions v n defined by 

v n (x) = tn (^7|) expii() n (x), 

where 

f 0, if x> i, 

ipn(x) = I Qn(nx - 1), if - i < x < i, 
[ -2q n , if x < -i, 

with n > 1 and = -, r , and compute 

E(Vn) _ ^ + 2^ - th(^-)) - ^, as „ + oo, 



whereas 
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In order to prove Theorem 2, we need to analyse minimizing sequences for (Vp). More 
generally, for given < p < we consider a sequence (u n ) ne ^ verifying 

p n = p(u n ) -> p, and E(u n ) -> E min (p), as n ->■ +00. (2.27) 

Minimizing sequences for E m [ n (p) are a special example of sequences satisfying (2.27). We have 

Theorem 3. Let < p < ^ be given, and let (ii n )neN be a sequence of maps satisfying (2.27). 
Then, there exist a subsequence (u a r n )) n ^, a sequence of points (a n ) ng N, and a real number 8 
such that 

u a (n){ • +o<T(n)) -> expi6» C ( P )(-), as n -> +00, 
uniformly on any compact subset o/R. Moreover, 

1 - K(n)(' +M«))| 2 ^ 1 - I f c(p ) (-) I 2 in L 2 (R), asm +00, 

and 

%„,)(■ +a CT ( n )) expi6> o' c(p) (-) in L 2 (R), as n -> +00. 

The proof of Theorem 3 relies on several mostly uncorrelated observations, which we next 
present as separate lemmas. The first elementary observation emphasizes the role of the sonic 
speed y/2. 

Lemma 4. Let g and tp be real-valued, smooth functions on some interval of M, such that g is 
positive. Set v = gexpiip. Then, we have the pointwise bound 

(1 - g 2 )ip' < —e(v). 

Q 

Proof. Notice that we have by (2.18), 

<v) = \({g') 2 + g\v'?) +\{l~Q 2 ) 2 > \(e 2 W? + 5(1 - Q 2 ? 

The conclusion follows from the inequality \ab\ < \(a 2 + b 2 ) applied to a = -^(1 — g 2 ) and 
b = gtp'. □ 

As a consequence, we have 
Corollary 1. Let < p < 5, and v S X p . Then, 

E(v) 



inf \v(x)\ < 

In particular, if 5{v) = 1 — ^jS^ > 0, then, given any < 6 < S(v), there exists some point 
xs G K such that 

1- \v(x s )\ > S- 

Proof. Set £0 = influx)!. It follows from Lemma 4 that we have the pointwise bound 



(i - e 2 y 

and the conclusion follows by integration. □ 



< — e(v) 
QO 
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The next observation shows that the size of regions where \v\ is away from 1 can be bounded 
in terms of the energy. 

Lemma 5. Let E > and < 8q < 1 be given. There exists an integer £q = £o(E, So), depending 
only on E and Sq, such that the following property holds: given any map v G (M) satisfying 
E(v) < E, either 

|l - \v{x)\ \ < S , Vx G R, 
or there exists I points X\, x-i, ■ ■ ■, and X£ satisfying i < Iq, 

|l - \v{xi)\\ >5 , VI < i < £, 

and 

1 - \v(x)\\ <5 , Vx £ R \ U [Xi - l,Xi + ll. 

i=l 

Proof. Set 

A = {zeR, s.t. |l - |^(^)| | > 5 }, 
and assume that A is not empty. Considering the covering R = .U^ii where Jj = [i — ^, i + 
we claim that, if lj fl A ^ 0, then 

e(«) > Mo, (2.28) 

where Ij = [i — 1, i + 1], and /xo is some positive constant. To prove the claim, we first notice 
that 

\v(x) - v(y)\ < \\v'\\ L 2 m \x - y|3 < y/2E*\x - y\* , (2.29) 
for any (x, y) G R 2 . Therefore, if z £ A, then, 

|l - \v(y)\\ > —, yy G \z - r,z + r], 

<5 2 i 
where r = Choosing ro = min{r, g}j we are led to 

e{v)> \ / (1-M) 2 >^^. 



4 , 

z— ro J z— ro 



In particular, if z G ij Pi .4 for some i G N, then [z — ro, 2 + ro] C lj, and claim (2.28) follows. 
To conclude the proof, we notice that 



so that, in view of (2.28), 



tHo < 2E, 



where I = Card{z G N, s.t. Ij n .A / 0}. The conclusion follows setting £ = — = 12 ff , and 



/'" 



choosing some point Xj G fl .4, for any i G N such that Ij fl .4 7^ (relabelling if necessary, the 
points Xi). □ 

We will also need the following construction. 
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Lemma 6. Let < |q| < ^ and < p < \. There exists some number i > 1, and a map 
w = \w\ exp iip G -ff 1 ([0, £]), such that 



w (0) = w (£) : < 1 1 — |io(0)|| < /i, 



and 



#M < 14|q|. 



(2.30) 
(2.31) 

(2.32) 



Proof. Consider the functions f% and tp\ denned on the interval [0, 2] by 



r l-i rl I 

/l(s) = s on |0, -j , /i(s) = 1 - s on lj , and /i(s) = on 



1,2 



and 



V'i(s) = s on [0, 1], and i/ji(s) = 2 — s on [0, 1]. 
For a given positive number A > 0, we consider the functions defined on [0, 2A] by 

/a(s) = -^/(j), and tp x (s) = V'Q), 
so that |/ A | < \i>' x \ = i / A (0) = / A (2A) = 0, Va(0) = </>a(2A) = 0, and 



2 A 



l 

4A' 



2A 



4 : 



2A 



fx 



1 

12A : 



2A 



1 

A^ : 



(/a) 2 = and 



2A 



W>a) 2 = | 



We then choose A = ^-r, so that, j < introduce a new parameter 5 > to be determined 
later, and consider the function 



Px,5 



so that 1 — p\ g = f\ + 5. It follows from our choice of parameter A that 



2A 

fx^'x = o 



2A x 

(/a + <M = - 



2 A 



(i - pI^a- 



(2.34) 



We finally choose I = 2A and 



•«' 



px,sexpiip\, if q > 0, 
Px,&ew(-i^x), if q < 0. 



Condition (2.31) is fullfilled with this choice of w in view of (2.34). Moreover, by construction, 
w(0) = w(£) = px.s = v 7 ! — S, so that conditions (2.30) are satisfied for any 5 < p? . We finally 
compute 



E(w) 



2A 



(fx) 2 



8(1-5- fx 



+ {1-5- f. 



\ m 2 , fi , sfx . 

7 2 4 2 4 



so that, since 



it follows from (2.33) that 



°^ + ^2T + ^ 2 n4 



E(w) < 



2A /(fl\2 



(fx) 2 , W , fx , 5/a , 



+ 



+ — + + 
2 4 2 4 



1 1 1 5 5 2 X 
-4A3 + A + 48A + 8 + ~2~' 



Inequality (2.32) follows choosing 5 = min{/i 2 , j}. 



□ 
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We are now in position to undertake the study of sequences satisfying (2.27). We set 

E(u n ) 



S(u n ) = 1 



V2p(u n ) ' 



so that, by (2.27), 5 n — > 5 P = 1 — , as n — > +oo. Moreover, it follows from the properties 

of the curve p i— > £(p) (see Subsection 2.2) that 

5 P > 0, (2.35) 

for any < p < |, so that 

S n > |, (2.36) 

for n sufficiently large. We first have 

Lemma 7. Let < p < f, and /e£ (?x n ) ng N be a sequence of maps satisfying (2.27). Then, there 
exists a subsequence (u a ^) ne ^, some numbers 6 and x, and a solution v c to (TWc) such that 

u a(n) expi# d c (- + x) weakly in i? 1 ([— A, A]), as n — > +oo, 
for any A > 0. In particular, we have 

[ A e(t> c ) < liminf f /" e(u n )) , and (1 - £)<f/ c = lim f [ (1 - , ( 2 . 3 7) 

J -A ™->+oo\J_ A J J_ A n^+oo\J_ A J 

where we have written u n = g n expi(p n , and X> c = g c expi<p c . 

Proof. Since (E(u n )) n ^ is bounded by assumption (2.27), it follows from standard compactness 
results that there exists a subsequence (u a ^) n ^, and a map u G flj^R) such that 

u a(n) u 111 H 1 ^— ^4)^4])) as n — > +oo, 

for any ^4 > 0. Assertions (2.37) follow using Rellich's compactness theorem, and the lower 
semi-continuity of E. Hence, it remains to prove that the limiting map u solves (TWc) on 
(—A, A). For that purpose, we consider a smooth map £, with compact support in (—A, A), 
such that 

(m,£')=0. (2.38) 



We claim that, for any t sufficiently small, 

f e{u a{n) +tt)> I' e(u a(n) ) + 0(t 2 ) + o(l). (2.39) 

J — A J— A n— >oo 

To establish the claim, we first observe that, in view of (2.27), we may assume without loss of 
generality that E(u n ) < ^p, so that u n has no zero. Expandind p{u n + 1£), we obtain 

p(u n + tf ) = p(u n ) + t / (iu n ,£'} + 0{t 2 ) 

= Pn + 0(t 2 ) =p + 0(t 2 ) + o(l), 

n^oo 

so that, setting q n t = p — p(tt n + t£)> we are to 

q n = 0(t 2 )+o(l). 
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We next construct a comparison map v n t for E mill (p) applying several modifications to the map 
Un + t£. For that purpose, we invoke Lemma 6 with q = q n j, and \i = fi n ,t = inf{-|, 
where u n> t = sup{|l — [—A, A]}. This yields a positive number £ n ^ > 1, and a map 

w n ,t = |^n| expi^ n , defined on [0,£ n (t)] such that 

w n (Q) = w n (£ n>t ), and |l - |w n , jf (0)|| < fj, n>t , 

and such that 



1 

In,* = 

and 



In,* = | / (1 - |^„,t|) 2 <, 







< 14|q„, t | = 0{t 2 ) + oil) . (2.40) 

n—*oo 

In view of the mean value theorem, there exists some point x n in [A, +oo) such that |w n (a; n )| = 
|^n,t(0)|. Multiplying possibly w n ^t by a constant complex number of modulus one, we may 
therefore assume, without loss of generality, that u n {x n ) = w n (Q). We define the comparison 
map v n j as follows 

V n ,t( X ) = u n{x) + t£(x), Vx < X n , 

v n>t (x) = w n (x - x n ), \/x n <x<x n + £ n>t , (2.41) 

V n ,t(. x ) = u n{x - i n ,t) + t£(x - £ njt ), \/x>X n + £ n>t . 

We verify that 

E{v nt t) = E{u n + tt) + E(w ntt ), and p(v n>t ) = p(u„ + t£) + q n , t = p, (2.42) 
so that v n ^ is a comparison map for E m \ n , and therefore 

E(v nt t) > E min (p). (2.43) 
On the other hand, we have in view of assumption (2.27), 

E(u n ) = £ min (p) + o(l), (2.44) 
whereas, since £ has compact support in (—A, A), 

E{u n + tt)-E{u n ) = J (e(u n +t£)-e{u n j). (2.45) 

Combining (2.45) with (2.44), (2.43) and (2.40), we establish claim (2.39). 
To complete the proof of Lemma 7, we expand the integral in (2.39) so that 



t 



[ (<£' + e«n(i - Kl 2 )) > o(t 2 ) + 0(1) . 

J— A n— >oo 



We then let n tend to +oo. This yields, in view of the compact embedding of H l {[—A,A\) in 
C°([-A,A]), 



t J + £«(1-M 2 )) >0(t 2 ). 



Letting t tend to + and , we deduce 

-A 



[ u'£' + £u(l - \u\ 2 ) = 0. 
J-A 



Since £ is any arbitrary function with compact support verifying (2.38), this shows that there 
exists some constant c such that u solves (TWc), and hence is of the form u = expiO t> c (- + x), 
for some c > 0. □ 
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Remark 6. At this stage, we might only consider maps o c with positive speed c, using the 
convention that, for any c < 0, c (x) = D_ c (— x). 

It might happen that the limit map provided by Lemma 7 is a constant one. To capture the 
possible losses at infinity, we need to implement a concentration-compactness argument. The 
first step in this argument is to invoke Lemma 5, and Corollary 1 to assert 

Proposition 2. Let < p < \ be given, and let (it n )neN be a sequence of maps satisfying (2.27). 
There exists an integer £ p , depending only on p, such that there exists £ n points x\, x^, . ■ ., and 
x™ satisfying £ n < £ p , and 

|l-Mx?)|| > j, VI (2.46) 

and 

|l-K(x)|| VxeR\ U \x?-l,x? + l). (2.47) 

4 i=l 

Passing possibly to a subsequence, we may assume that the number £ n does not depend on 
n, and set i = l n . A standard compactness argument shows, that passing again possibly to a 
further subsequence, and relabelling possibly the points x™, we may find some integer 1 < I < £, 
and some number R > such that 

\xi - x]\ -> +oo, asm +oo, Vl<i^j<£, (2.48) 

and 

x? e U B(xf, R), W<i<£. (2.49) 
Going back to Proposition 2, we deduce 

5 I 

|l - \un(x)\ \ < 4, E M\ U B(x?,R + 1), (2.50) 

4 i=l 

so that, invoking Lemma 4, we have on M. \ U B(xf, R + 1), 



i=l 



e(ii n ) 



We are now in position to provide the proof to Theorem 3. 

Proof of Theorem 3. We divide the proof into several steps. 
Step 1. For any 1 < i < £, there exists numbers Ci ^ 0, x~i and B{ such that 

u n {- + xf) -± exp(i6i)t) Ci (- + Xi). 

Applying Lemma 7 to the sequence u n (- + x™) ng N yields the existence of the limiting solution 
to (TWc), exp(i6i)o Ci (-+x~i). It remains to prove that the function t) Cj is not a constant function. 
This is a consequence of the fact that 

\u n (x2)\ < 1 - ^ 

so that, since by compact embedding, we have uniform convergence on compact sets, we obtain 



|t>c,(fi)|<l--f. 



Hence, t> c . is not a constant map. 
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Step 2. Given any number [i > 0, there exists a number > 0, and £ N, such that, if 
n> no, then 



1=1 



e(u n ) > ^£(o c ,) 



anc 



i — 1 



(el - iK - J> 



< 



//.. 



where pi = p(o c 



To prove Step 2, we choose A > R + 1 so that, for any 1 < i < I, we have 

'■ A ii 

e(v Ci ) > E{v Ci ) - -~, 
-A 2£ 

and 



(\v c f - l)tp' c - p, 



< 



The conclusion follows from the convergences stated in (2.37). 
Step 3. We have 



U B(x",A) 



< 



1 



V2(l - "f) MvB(xf,A) 



e(u n ). 



To establish this inequality, it is sufficient to integrate (2.51). 

Passing possibly to a further subsequence, we may assume that there exist some numbers p^ 



and such that 



7: / (el ~ l)^n -> Pm> and / i e ( u n) -> EfM, as n -> +oo. 



U ^(x™,^) 



Going back to Step 2, and letting ?i — > +oo, we are led to the estimates 



i=l 



p - Y^Pi -Pv < A*, and 8(p) > Y^ £ (Pi) + E ^~^ 



i=l 



whereas Step 3 yields 

V2(l-^) Pfl <E,. 

Letting [i — > 0, we may assume that for some subsequence (/x m ) m eN tending to 0, we have 

P^m ~> P> and E nm -> E > as m -> +oo. 
Our previous inequalities then yield 



P = j> 



i=l 



(2.52) 



i=l 



with 



^(i-a)><*. 



(2.53) 
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Step 4. We have E = p = 0, and £ = 1. 

This statement is a consequence of properties of the curve p i— > £(p), and the definition of 
Assume first that p < 0. Then, it follows from (2.52) that 

I 

P (2-54) 

i=l 

and 

£(p) > J>(N). (2-55) 

i=l 

On the other hand, as observed before, the function £ is concave on M + , non-decreasing, and 
£ (0) = 0, so that £ is subadditive, and hence in view of (2.54), 

I 

£(p) < J2 £ (1*1). 

i=l 

so that (2.55) becomes an equality. However, since £ is strictly concave, this is possible if and 
only if I = 1, and (2.54) is an equality. Step 4 therefore follows in the case considered. 

It remains to handle the case p > 0. In view of the definition on 5p, we have 

^(P) _ . /oti X \ ^ . /n f i S ? 



P 

so that by concavity, we have for any < s < p, 



V2(1-5 P ) < ^(l-^), (2.56) 



£{s) > s^- = sV2(l - 5 P ) . (2.57) 



Setting p = p - p = ^2£(pi), we obtain using (2.52), (2.53), (2.56) and (2.57), 



i=l 



> £(p) - V2p{l -*±)> £(p) - E 



(2.58) 



> 



£*(i* 



i=l 



We then conclude following the first case, replacing p by p, and using the subadditivity of the 
function £. 

Step 5. Proof of Theorem 2 completed. 

At this stage, we have shown that there exist some 6 = 9±, and some solution r> c = D C1 such 
that 

u n (- + On) expiO tt ci , 

weakly in iif 1 (lR), as n — ► +oo, where a n = x™ — x~i, and that we have the convergence of the 
energy E{u n ) — ► £(b Ci ). In particular, 

■u' n (- + a n ) — > expz# \}' C1 , 

strongly in L 2 (R) as n — ► +oo, and the conclusion follows by compact embedding. □ 



19 



2.4 Orbital stability 

We recall first the classical notion of orbital stability (see, for instance, [5]). A solution D c to 
(TWc) is said to be orbitally stable in a metric space -X'(M), if and only if given any e > 0, there 
exists some 5 > such that for any solution ^ to (GP) in X(R), if 

djc(*(-,0),d c ) < 5, 

then 

sup I inf dx(^(-,t),expi8 C (- — a)) ) < e, 

tm \ (a,e)m 2 ) 

As a preliminary step, this definition requires of course to prove first that the Cauchy problem 
for (GP) is globally well-posed in X(R) (see [30]). A natural choice is the energy space 

X 1 (R) = {w £ L°°(R), s.t. w' G L 2 (R), 1 - \w\ 2 G L 2 (R)}. (2.59) 

Given any vq G ^(R), Zhidkov [67] (see also [30]) established that (GP) has a global solution 
^ with initial datum vq. Moreover, the Ginzburg-Landau energy is conserved, 

£?(*(-, t)) = E(yo), Vt G R. (2.60) 

Notice that E is a continuous function on ^(R). If we assume moreover that 

E(v ) < ^, 

using (2.60), we are led to 

E(*(;t)) < Vi G R, 

so that, by Remark 1, we may write \t(-, t) = g(-, t) exp i(p(-, t) for any ( £ 1. We may therefore 
define the scalar momentum p(^>(-, t)) using (2.22). It is shown in [50] that the scalar momentum 
p is then a conserved quantity, i.e. 

p(* (•,*)) =p(v ), Vt G R. (2.61) 
Given any A > 0, we consider on A' 1 (M) the distance x 1 defined by 

dA,X^{vi,V 2 ) = \\vi - «2 1| !,«»([- A, A]) + IK - ^2lli 2 (K) + IIKI ~ l D 2| WtfQg) ■ ( 2 - 62 ) 

Following ideas from Grillakis-Shatah-Strauss [39], Zhiwu Lin [50] proved the orbital stability 
of the solutions D c to (TWc), for any < c < y/2, when the perturbations w are taken in the 
space X 1 (R). 

Using our study of the minimization problem (Vp), we would like to give a short proof of the 
result of Zhiwu Lin [50] which follows the classical compactness argument of [17]. We first recall 

Theorem 4 ([50]). For vq G ^ 1 (M), consider the global in time solution ^ having initial datum 
vq. Let < c < y/2 be given. For any numbers e > and A > 0, there exists some positive 
number 5, such that, if 

^i(»o,t) c ) <6, (2.63) 

then, we have 



sup inf_ d Ax i^(-,t),expi9 B c (- - a)) ] < e. (2.64) 
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Proof. We assume by contradiction that D c is not orbitally stable for the distance d A x 1 ■ ^ Q 
this case, we may find a positive number eo, sequences of numbers (5 n )«eN and (tn)nGN, and a 
sequence of functions (t>o)neN such that 5 n — > 0, as n — > +00, 



d 



(«o,t> c ) < S n , 



and 



inf d A x i{ty n (-,t n ),expi9 C (- - a n )) > s , 



(2.65) 



(2.66) 



where \E ,Tl denotes the solution to (GP) with initial datum Vq. Denoting w n = *$> n (-,t n ), the 
conservation of the Ginzburg-Landau energy and of the scalar momentum implies that 



E{w n ) = E(v%), and p(w n ) = p(v%), 
for any n £ N. Moreover, the energy is continuous on X 1 (R), so that, invoking (2.65), 

E(vjj) -> E(x> c ) = E min (p), as n -> +00, 
where p = p(t) c ) is given by formula (2.24). Therefore, by (2.67), 

E{w n ) -> E min (p), asm +00. 
In view of (2.22), the scalar momentum may be written as 



(2.67) 



p(v) 



~-{iv,v'), 



for any function v 6 A' 1 (R). By (2.65), we have 



), as n — > +00. 



(2.68) 



(2.69) 



(2.70) 



On the other hand, in view of (2.4), the modulus of o c has a minimum value equal to on M 
which is achieved at the origin. Since (\vq |) n eN uniformly converges to |o c | on M as n — > +00 by 
(2.65), we may assume that 

inf\vS(x)\ > 



2^2' 



for n sufficiently large, so that 



,,n|2 



,n|2 



4 



,n|2 



By (2.65), it follows that 



Kl 2 - 1 

\n,n\2 



|d c | 2 — 1 2 . 

— — — in L (K), as n — > +00, 



so that by (2.69), and (2.70), 



p(«o ) ~> -P(°c) = P, as n ^ +00. 



(2.71) 



By (2.67), we are led to 

p(w n ) — >■ P, as n — ► +00, 

so that, by (2.68), the sequence (w; n ) n£ N verifies assumptions (2.27). Hence, by Theorem 3, 
there exist some points (a n ) ng N, and some real number 9 such that, up to some subsequence, 



w 



n{ ■ +a n ) ->• expz<9 C ( P )(-), as n ^ +00, 



(2.72) 
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uniformly on any compact set of K, and 

w 'n{ ' +°n) — ► expi# t)^(-) in L 2 (R), as n — > +00. (2.73) 
In view of (2.68) and (2.73), we have 

/ (1 - \w n \ 2 ) 2 ->■ / (1 - |o c | 2 ) 2 , as n > +00, 
Jr Jr 

whereas, by (2.68) and (2.72), up to some subsequence, we have 

1 - \w n ( • +a n ) | 2 — 1 — |t) c ( p )(-)| 2 in L 2 (R), as n — > +00. 

This yields 

1 — |u>n( ' +«n)| — ► 1 — |&c(j>)Q| in L 2 (M), as n ^ +00. 
Hence, by (2.72) and (2.73), 

Wn( ■ +fln) — ► expi# tJ c ( p )(-) in A^M), as n — > +00, 

which gives a contradiction with (2.66), and completes the proof of Theorem 4. □ 

Remark 7. Di Menza and Gallo [25] proved the linear stability of Do submitted to small per- 
turbations in {d } + H 1 ^) (see also [30]). 

2.5 Relating (TWc) to the Korteweg-de Vries equation 

Travelling wave solutions to (TWc) are related to the soliton of the Korteweg-de Vries equation 
as follows. Set e = \/2 — c 2 , and consider the scaled function 

N e (x) = ^ V J^\ (2.74) 



e 2 -\e 



where r\ c = \ — |o c | 2 . Invoking (2.16), we are led to 

1 



N s (x) = N(x) 



2ch 2 (f) 



A remarkable property of N is that it represents the classical soliton to the Korteweg-de- Vries 
equation 

d t w + dlw + 6wd x w = 0. (KdV) 
Concerning the phase (p c of D c , we consider the scale change 

Q e (x) = ^<p c (-\ (2.75) 



e \e 

so that we similarly obtain from (2.17), 



It is also of interest to compute the corresponding energies. Using scale changes (2.74) and 
(2.75), we are led to 

„, , e 5 f (N') 2 e 3 / e 2 \ f N 2 e 3 f , 



22 



and 

so that the quantity \/2p(t) c ) — £7(t> c ) is governed, at the limit e — ► 0, by 

e 5 

V2p(* c ) - E(X> C ) ~ -—E KdV {N), 
where Exdv{^) is the energy of the soliton to (KdV), namely 

iWAO = i / {N'f - / 

1 Jr Jm. 

3 Qualitative properties of travelling waves in higher dimen- 
sions 

We next turn to finite energy travelling waves in higher dimensions, and describe a number of 
qualitative properties, which have been rigorously established so far. Many of the results in this 
section were already guessed in the quoted seminal papers of Jones, Putterman and Roberts 
[46, 45]. Some mathematical proofs turn out to be quite different from the physical intuition. 

3.1 Range of speeds and energies 

The following results are proved in [34, 36]. They describe the possible spectrum of speeds c, 
excluding in particular the possibility of non-constant supersonic travelling waves. 

Theorem 5 ([34, 36]). i) Let N > 2. Any travelling wave of finite energy and of supersonic 
speed c > y/2 is constant. 

ii) Let N = 2. Any travelling wave of finite energy and of sonic speed c = v2 is constant. 

We briefly sketch the proof of Theorem 5 in the supersonic case. The proof is similar in the 
sonic case. 

Proof. Let v be a finite energy solution to (TWc) of speed c > y/2. In order to prove that v 
is a constant function, we compute some integral identities (relating the energy and the scalar 
momentum), so that the energy of v is necessarily equal to 0, i.e. v is a constant function 
of modulus one. The first identities are the so-called Pohozaev's identities [58], obtained by 
multiplying (TWc) by the test function (x,Vv(x)} , and integrating by parts. 

Lemma 8. Let N > 1 and c > 0, and consider a finite energy solution v to (TWc). Then, 

E{v) = [ \d lV \ 2 , (3.1) 

and for any 2 < j < N , 

E(v) = I \djv\ 2 + cp{v). (3.2) 

Identities (3.1) and (3.2) are not sufficient to ensure that the energy of v is equal to 1 , so 
that another integral relation is required. 

1 Notice that (3.1) and (3.2) are sufficient to ensure that any travelling wave of finite energy and of speed c = 
is constant (see [46, 11]). 
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Lemma 9. Let N > 2 and c > \/2, and consider a finite energy solution v to (TWc). Then, 

/ kji (|V^ + (1-H 2 ) 2 )=2c(i-|> W . (3.3) 

Identity (3.3) is obtained using the singularities of the kernels associated to (TWc) in case 
c > y/2. Indeed, assuming that v does not vanish on M. N (in order to simplify the proof), v may 
be written as 

v = gexp iip, 

where g = \v\, and cp are smooth, real- valued functions. By (TWc), the function n = 1 — g 2 is 
solution to 

A 2 n - 2An + <?d\ x r\ = —AF + 2cdidiv(G), (3.4) 

where 

F = 2\Vv\ 2 + 2rf - 2crjd 1 ip, and G = -ifJip. (3.5) 

Notice in particular that the nonlinearities F and G are related to the density of energy and 
momentum (equal to — G in view of (2.22)). Taking the Fourier transform of (3.4), we are led 
to 

N 

m = K (OH)-2cJ2 K i(OG(O, (3.6) 

J'=l 

where K Q (£) = |g|4 +2 |||2_ c 2g2 , and Kj(£) = |g|4 +2 ^p_ c 2g2 , for any 1 < j < N. The singular 

nature of the kernels Kq and Kj now gives some relation between F(0) and G(0). This in turn 
yields a relation between the energy and the scalar momentum. Indeed, since v is of finite 
energy, the function rf in the left-hand side of (3.6) belongs to L 2 (R N ), whereas the functions 
F and G in the right-hand side are continuous. When < c < y/2, Kq and Kj are sufficiently 
smooth, so that (3.6) can hold without additional assumption. In contrast, when c > \/2, Kq 
and Kj have singularities on the set 

r = {£e]RM£l 4 + 2|£| 2 -c 2 ^ = o}. 

Since (3.6) holds, this leads to a relation between and G(£) on the set T. Taking some 

limit ^ — > 0, this yields identity (3.3). Combining (3.1), (3.2) and (3.3), it may be shown that 
the energy of v is equal to 0, so that v is a constant of modulus one. This completes the proof 
of Theorem 5. □ 



In the three-dimensional subsonic case, small energy solutions are also excluded in view of 
the following theorem (see [9]). 

Theorem 6 ([9]). Let N = 3 and < c < s/2. There exists some positive universal constant £q 
such that any non-constant finite energy solution to (TWc) on M 3 satisfies 

E{v) > S . 

Theorem 6 improves an earlier result by Tarquini (see [65]), which states that any solutions 
of sufficiently small energy, with respect to their speed, are excluded in any dimension N > 2. 

Theorem 7 ([65]). Let N > 2 and < c < y/2. There exists some positive constant S(c,N), 
depending only on c and N, such that any non-constant finite energy solution to (TWc) on M 3 
satisfies 

E(v) > £{c,N). 

Moreover, 

£{c,N) ->0, as c ->■ y/2. 
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Theorem 6 leaves some hope for a complete scattering theory for solutions with small energy 
in the three-dimensional case. Such a theory has been established in any dimension N > 4 by 
Gustafson, Nakanishi and Tsai (see [41], or [56] in the present volume). In dimension three, 
they were able to establish the existence of global dispersive solutions to (GP) (see [42] or [56] 
in the present volume). On the other hand, two-dimensional travelling waves of small energy 
are known to exist (see Theorem 13), excluding scattering in the energy space. 

For sake of completeness, we briefly sketch the proof of Theorem 6. 

Proof. The proof relies on equation (3.6), using the next elementary observation. If some quan- 
tity E satisfies 

E < KE 2 , (3.7) 
where K is some positive constant, then E is either equal to 0, or 

E >I. 

~ K 

Taking the L 2 -norm of (3.6), we are led to an inequality of the form (3.7), where the energy 
E(y) plays the role of the quantity E, whereas K is equal to some linear combination of the 
L 2 -integrals of the kernels Kq and Kj . 

Indeed, first notice that the nonlinearities F and G in the right-hand side of (3.6) are (al- 
most) quadratic quantities with respect to rj and Vip, related to the densities of energy and 
of momentum. In particular, their L 1 -norms are bounded by the energy up to some positive 
universal constant K, so that 

\F\ + \G\ < KE{v). 
Taking the L 2 -norm of (3.6), we are led to 

II^/IIl2(rjv) = ||??||l2(rjv) < ke(v)\\Kq\\ L 2 



On the other hand, it is proved in [9] that any finite energy solution v to (TWc), with the 
additional assumption \v\ > \ (which holds in case E(v) is sufficiently small (see [9])), satisfies 

E(v) < 702117/1112 

so that, by the above elementary observation, 

E(v) > — . (3.8) 

7c 2 K^\\K \\l 2(m) 

A direct computation now gives 

ll-Kolliapis) = ^ arcsin (-^=) , 

so that 

E( s > 1 > £ = — 

U ~ 7vr 2 i^ 2 c arcsin (-^) ~ ° 7vr 3 K 2 ' 

This completes the proof of Theorem 6. □ 

Remark 8. The proof above cannot be performed in the two-dimensional case. Indeed, a direct 
computation gives 

\\ k o\\l 2 



V2(2 - c 2 ) 
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so that (3.8) now becomes 



E{v) * £{c) = rS^-> (3 ' 9) 

where £ (c) — > 0, as c — > \/2- In particular, (3.9) does not prevent the existence of solutions with 
energy as small as possible (see Theorem 13). 

3.2 Regularity and decay at infinity 

The following results were proved in [26] and [9] (see also [11, 35, 65, 64]). They describe the 
regularity of subsonic travelling waves (in particular their real-analyticity) . 

Theorem 8 ([26, 9]). Let N > 2 and < c < a/2, and consider a finite energy solution v to 
(TWc). Then, v is a real- analytic, bounded function on M. N such that 



M<yi + j- (3-10) 
For sake of completeness, we briefly sketch the proof of Theorem 8. 

Proof. The smoothness of v follows from a standard bootstrap argument using the finiteness of 
the energy, and the elliptic nature of (TWc). Bound (3.10) essentially results from the maximum 
principle for the function \v\ 2 , which verifies 

|2 , ol„.|2/i , c |„.|2\ _ olV7„.|2 o„/„fl „. „.\ , 6 L.I2 



A\v\ 2 + 2\v\ z \l + — - \v\ z j = 2\Vv\ 2 - 2c(idiv ,v) + — \v\ z > 0, 

by (TWc) (see [26, 65, 64] for more details). 

Real-analyticity is established following an argument of Bona and Li [13, 14] (see also [47, 
52, 53]). The idea is to prove the uniform convergence of Taylor's series of v, 

d a v(x) 
~~d~\ 



on a complex neighbourhood of any arbitrary point x 6 R . The required estimates for the 
derivatives are provided by (TWc), using standard L 9 -multiplier theory, Sobolev's embedding 
theorem, Gagliardo-Nirenberg's inequality, and the superlinear nature of the nonlinearities. 

Indeed, denoting v\ = Re(v) — 1 and v 2 = \m(v), equation (TWc) may be recast as 

dj k V! = H jjk * F 1 (vi,v 2 ) - icHi jjik * F 2 (v 1 ,v 2 ), (3.11) 
dj k v 2 = icH ltjsk * Fx(vi,v 2 ) + K j)k * F 2 (vi,v 2 ), (3.12) 

where the nonlinearities F\ and F 2 are defined from C 2 to C by 

Fi(z\, z 2 ) = 3zj + z\ + z\ + ziz 2 , and F 2 (zi, z 2 ) = 2z\z 2 + z\z 2 + z\ , 

and the kernels Hj tk , Hij^ and Kj jk are given by 

fTtn 7T~(f\ - htek fT(ft= ^ + \^\ 2 ) 

*jMV |^|4 + 2 |^|2 _ c 2£2 > |£|4 + 2 |f|2 _ < A jMU |£|4 + 2 |£|2 _ > 

for any 1 < j,k < N. By a result from Lizorkin [51], the kernels Hjk, Hi,j,k an d Kjk are 
L^-multipliers for any 1 < q < +oo. Coupled to (3.11) and (3.12), this provides L 9 -estimates of 
any derivatives of v\ and v 2 . 
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Lemma 10. Let 1 < j,k < N, a G N w and 1 < q < +oo. There exists some positive number 
K\(q), possibly depending on q, but not on a, such that 

H^^iHi^r") + W^^k^hi^ < ^(^(ll^Fi^i,^)^^^) + \\d a F 2 (v 1 ,v 2 )\\ Lq{ 

(3.13) 

The second step is to transform L^-estimates (3.13) into uniform ones. This follows from 
Sobolev's embedding theorem and Gagliardo-Nirenberg's inequality. 

Lemma 11. Let 1 < j < N, a € and y < q < +oo. There exist some positive numbers 
K 2 (q) and K 3 (q), possibly depending on q, but not on a, such that 

\\d a v 1 \\ Lac{R N ) + \\d a v 2 \\ L o C(R N ) <K 2 (q)F q (a), 

(3.14) 

\\d a djVi\\Lci(R") + H<9 Q djW2||z/j(iRiV) <K 3 (q)F q (a), 

where we have set 

F q (a) = m&x^Wd^F^vuv^W^^N) + \\d^F 2 {v u v 2 )\\ Lt ^N^j . 

In view of (3.14), the convergence of Taylor's series T VljX (z) and T V2iX (z) reduces to the 
convergence of the series 



al 



\Z - X 



for z sufficiently close to xq, and q suitably chosen. Combining the superlinear nature of F± and 
F 2 with estimates (3.13) and (3.14), an inductive argument based on Abel's identity as in [52], 
yields 

Lemma 12. Let a S and y < q < +oo. There exists some positive number K^q), possibly 
depending on q, but not on a, such that 

F q (a) < K 4 (q)Ma", (3.15) 

where we have set a = (maxjai — 1, 0}, . . . , maxjctAr — 1, 0}). 

In view of (3.15), choosing q = N, the series Sg )Xo (z) is convergent for any z such that 
\z — xo\ < K j N \ ■ Taylor's series T Vl ^ Xo (z) and T V2tXQ (z) converge the same way, so that v is a 
real-analytic function on M. N . □ 

Jones, Putterman and Roberts [46, 45] investigated the decay properties of subsonic travelling 
waves in dimensions two and three. They computed some formal asymptotics of axisymmetric 
solutions to (TWc) in [46, 45]. Their formal derivation was confirmed, and somewhat extended, 
in [33, 35, 37, 38]. 

Theorem 9 ([33, 35, 37, 38]). Let N > 2 and < c < and consider a finite energy solution 
v to (TWc). There exist a complex number Aqo, such that \Xoo\ = 1, and a smooth, real-valued 
function : S^ -1 — > M., such that 



Ixl^- 1 



(v{x) - Aoo) - iAoofoo -» in L 00 ^ -1 ), as \x\ -» +oo. (3.16) 



Moreover, there exist some real constants a, f5 2 , . . ., and Pn such that the function Voc is equal 
to 

N 

voc(o-) = a ^ + w , Va e S N ~\ (3.17) 

3=2 f 1 .^ + c_a 1 \ 
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The constants a and (3j are given by 

Pj = -^ 1 [ 1 -t) 2 p ^ (3 - 19) 

where Pj(v) = \ J Rn (idjV,v — 1). The constants /3j are equal to 0, when v is axisymmetric 
around axis x±. 

Notice that a finite energy solution v tends to some constant Aqo at infinity. Since multipli- 
cation by a constant of modulus 1 keeps (TWc) invariant, we may, without loss of generality, 
assume that A ro = 1. Notice also that the decay of the function v — 1 is algebraic. More precisely, 
v(x) — 1 decays as , ,^ f _ 1 , as \x\ — > +oo. This yields the following corollary of Theorem 9. 

Corollary 2 ([35, 37]). Let N > 2 and < c < \/2, and consider a finite energy solution v to 
(TWc). 

i) The function v — 1 belongs to L q (R N ) for any q > j^rj, and its gradient X7v is in L q (R N ) for 
any q > 1. Moreover, any higher order derivative d a v belongs to L q (R N ) for any q > 1. 

ii) Assume N = 2, and v is non-constant and axisymmetric around axis X\. Then, the function 
v — 1 does not belong to L 2 (M. 2 ). In particular, it does not belong to i7 1 (]R 2 ). 

Corollary 2 has some significant consequences. It first leads to a rigorous definition of the 
momentum P{v) of a finite energy solution v to (TWc). Indeed, formula (1.2) makes sense, 
in any dimension, for ^ = v , since v — 1 and X7v respectively belong to L 4 (M Ar ) and L3(M. N ). 
Actually, it makes sense for any to G W(R N ) = {1} + V(R N ), where V(R N ) is defined by 

V(R N ) = {v:R N ^C, s.t. (Vv,Re(v)) G L 2 (R N ) 2 , lm(v) G L 4 (R N ) and VRe(^) G L*(R N )}. 

(3.20) 



Notice that 
where 



W(R N ) C £(R N ) 



S(R N ) = {v£ Hl c (R N ), s.t. E(v) < +oo}, 



denotes the energy space. This last observation plays an important role in the variational 
argument of [9] for the construction of finite energy solutions to (TWc) (see Section 4 below). 
Our choice of the variational space is indeed W(R N ) = {1} + V{R N ). 

Statement ii) of Corollary 2 has to be considered in connection with the stability problem 
of two-dimensional travelling waves v. This analysis requires to find some functional space, on 
which the Gross-Pitaevskii equation is known to be globally well-posed, which preserves both the 
energy and the momentum, and which contains the travelling waves. As a matter of fact, there 
are several functional spaces where the two-dimensional Gross-Pitaevskii equation is known to 
be globally well-posed: for instance in {1} + ff 1 (IR 2 ) (see [11, 30]), in the energy space £(R N ) 
(see [31, 30]), and also in the space {w} + if 1 (]R Ar ) (see [28, 30]), for any w in the energy space 
£(R N ). Taking w = v , the space {v} + i7 1 (]R 2 ) seems appropriate to address the question of 
stability near a travelling wave v. Two important advantages of the space {v} + i? 1 (M 2 ) are 
that v belongs to this space (in contrast with {1} + H 1 ^ 2 ), in view of Corollary 2), and that 
the momentum is well-defined (in contrast with the energy space). 

For sake of completeness, we briefly sketch the proof of Theorem 9. 
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Proof. The proof is reminiscent of a series of articles by Bona and Li [13], De Bouard and Saut 
[24], and Maris [52, 53]. It relies on the use of convolution equations, in particular on a careful 
analysis of the kernels they involve. Indeed, assuming that v does not vanish on R (in order 
to simplify the proof), v may be written as v = gexpiip, where g and ip are smooth, real- valued 
functions. By (TWc), the functions r\ = 1 — g 2 and V</9 are solutions to (3.4), and 

A<^ = !&w + div(G), (3.21) 

where F and G are the nonlinearities given by (3.5). Taking the Fourier transform of (3.4) and 
(3.21), equation (TWc) reduces to the convolution equations 

N 

rj =K * F + 2c^ K 3 * G v 

N N 



k=l k=l 

where Kq, Kj, Lj /. and Rj /. are the kernels given by 



djip =^Kj *F + c 2 ^2 L i,k *G k + ^2 R 3,k * G k, 



(3.22) 



Both the decay properties and the asymptotics of v follow from the decay properties and the 
asymptotics of the kernels given by (3.22). The first step is to derive them. In view of (3.22), 
this reduces to the analysis of Kq and Rj,k- 

The kernels Rj )k are directly related to Riesz's operators Rj, given by 



R j(0 = 

since they may be written as 



Rj,k — —Rj * Rk- 

The properties of Riesz's operators are well-known (see [63]), partly because of their explicit 
expression 

\x\ 



7T 2 



where p. v. stands for the principal value. There is a similar expression for the kernels Rj t k: 

p M - v ( 6 ltk \x\ 2 - Nx jXk \ 

Rj,k{X) - K p.V. -- 1 y T 2 , [A.ld) 

2-k 2 V \M / 
so that their asymptotics are straightforward. 

In view of its non-homogeneity, the analysis of Kq is more involved. However, we have 

Lemma 13 ([35]). Let N > 2 and N — 2 < a < N. There exists some positive constant K a 
such that 

\K (x)\<^, Vx£R n . 
In particular, Kq belongs to L q (M. N ) for any jj-^j < q < jj^. 
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Lemma 13 results from some integrability estimates of the derivatives of Kq. They provide 
the algebraic decay of Kq by the inverse Fourier transform formula. Notice that this formula, 
coupled with some stationary phase estimates, states that the asymptotic properties of Kq are 
mainly given by the behaviour of Kq close to the origin. In view of (3.22), this behaviour is 
given by 

,£■2 -i ^ 

k (o ~ (2 _ c2 | + 2|c±|2 = (v^&, v~2u) + - 2 £ ^ (v^a, ^) . 

Up to some scale changes, the asymptotics of Kq are therefore the same as the asymptotics of 
Rjj. This gives for instance, 

Lemma 14 ([37]). Assume N > 2 cmci consider a smooth function f such that 

\f(x)\ < Y^w, and \Vf(x)\ < 1 + | ^ 2jV+1 , Vx G R w , (3.24) 
where K is some positive constant. Then, 

T(f )c 2 (2 - c 2 )^ (2 - c 2 + (c 2 - 2A0<r 2 

2x/2tt^ (2 - c 2 + c 2 a 2 



R n Kq * f(Ro-) -> — — — ^ — ^ ( / / ) , as R +00. (3.25) 



The asymptotics of a travelling wave v given by formulae (3.16), (3.17), (3.18) and (3.19) are 
derived using formulae like (3.25) for the kernels Kq, Kj, Lj^ and Rj^ki an d t ne nonlinearities 
F and G. Notice in particular that the presence of the energy and the scalar momentum in 
formulae (3.18) and (3.19) is due to the presence of the integrals of F and G in formula (3.25) 
above (see the proof of Theorem 5 above). 

Notice also that the proof of (3.25) first requires to establish (3.24), that is some algebraic 
decay for the function / convoluted to Kq, and its gradient V/. This is the second step of 
the proof. As in earlier papers [13, 24, 52, 53], this relies on some inductive argument using 
the superlinear nature of the nonlinearities F and G, which are almost quadratic functions of 
the variables 77 and in view of (3.5). When the nonlinearities are superlinear, the algebraic 
decay of a solution is identical to the algebraic decay of the kernels of the convolution equations 
it satisfies. 

To get a feeling for this claim, let us consider the simplified model 

f = K*f, (3.26) 

for which we have 

Lemma 15 ([35]). Assume r > 1, and consider smooth solutions f and K to (3.26), such that 
f belongs to L r (R N ), and K is in L 1 (R Ar ), and satisfies 

\K(x)\ < ^ — . 

' V >] ~ 1 + \x\ a K 

If there exists some positive constant < oq < ax such that 

\m\ < 



\x\ a ° ' 
then, 

\m\ < 



\x\ a K 

where A denotes some positive, possibly different, constants. 
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Lemma 15 states that the algebraic decay of / is exactly the same as the algebraic decay of 
K, if some small decay may be first established for /. Its assumptions are quite restrictive, but 
may be extended to more involved equations with only additional technicalities. The proof is 
by induction. Indeed, by (3.26), 



\x\ a \f(x)\ <A[ I \x- y\ a \K{x - y)\\f(y)\ r dy + / \K(x - y)\\y\ a \f{y)\ r dy _ 

(3.27) 

<Al ||| • r-K"|| Lt x=( R jvJ|/|| L T-( R jv) + H^HLim^)!!! ' I r ^IL° 



Using the assumptions of Lemma 15, equation (3.27) reduces to 

|l ' \ a f\\ roofloN\ < ^4 + ^4||| - | r /HioofTsW), (3.28) 



if < a < olk- Equation (3.28) now links the algebraic decay with exponent a of / to its 
algebraic decay with exponent ^. In particular, if we know some algebraic decay with a small 
exponent Qo > 0, a bootstrap argument yields the algebraic decay of / for a = qao, a = q 2 ao, 
. . ., that is for any a € [oq, a K\- This provides the optimal decay of /, which is the decay of the 
kernel K. 

However, in order to perform the inductive argument of Lemma 15, and to get the decay of 
a travelling wave v , we must establish some small decay for v. Multiplying equations (3.4) and 
(3.21) by r] and ip, and integrating by parts on the complementary of a ball, we may establish 
the following integral decay. 

Lemma 16 ([10, 35]). Let N > 2 and < c < \/2, and consider a finite energy solution v to 
(TWc) . There exists some positive constant qq such that 



e(v)(x)dx < +oo. 



Lemma 16 gives the small decay required to use the inductive argument of Lemma 15. Cou- 
pled with the decay and integrability properties of the kernels (see Lemma 13 and formula 
(3.23)), this first gives the decay of a travelling wave v , then its asymptotics, using formulae like 
formula (3.25) of Lemma 14. This completes the proof of Theorem 9. □ 



4 The existence problem in higher dimensions 
4.1 The variational approach 

As mentioned in the introduction, the existence problem in dimensions two and three has been 
widely investigated in the physical literature: rigorous mathematical proofs have been provided 
more recently, using a variational approach. 

It is a classical observation that one may obtain travelling waves by minimizing the energy 
E keeping the momentum p fixed. For a given p > 0, we therefore consider the minimization 
problem 

E min (p) = mf{E(v),v e W(R N ),p(v) = p}, 

where the definition of the space ^(M^) is given by (3.20). In this approach, equation (TWc) is 
the Euler-Lagrange equation to this constrained minimization problem. The speed c appears as 
a Lagrange multiplier, and is therefore not fixed a priori. Instead of using minimization under 
constraint, an alternate approach is to introduce, for given c, the Lagrangian 

F c (v) = E(v) - cp(v), 
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whose critical points are solutions to (TWc) . Solutions may then be found applying a mountain- 
pass argument. As a matter of fact, both techniques have been used so far. The first existence 
results were based on asymptotic Ginzburg-Landau's theory: the solutions obtained in that 
context had small speeds c > and possessed vorticity, i.e zeroes with non-trivial topological 
degree. For instance, in dimension two, one may obtain, using the mountain-pass lemma, a 
branch of solutions parametrized by the speed c > 0. 

Theorem 10 ([11]). Assume N = 2. There exists cq > such that there exists a non-constant 
finite energy solution v to (TWc) for any < c < cq. Moreover, v exactly has two vortices with 
degree ±1 located at a distance ~ -, as c — > 0, whereas 

2ir 

E(v) ~ 27rln(c), and p(v) ~ — , as c — ► 0. 





x 2 ' 

a 

c 


5 + i 






G 


Xl 

3-! 



In dimension three, the constrained minimization approach provides the following theorem 
for large momentums. 

Theorem 11 ([10]). Assume N = 3. There exists poo > such that, for any p > p^, there 
exists a solution u p to (TWc), with c = c(u p ), verifying p(u p ) = p, 

E(u p ) ~ TTv^h^p), and c(u p ) ~ as p ^ +oo. 

Moreover, v presents a vortex ring whose diameter is ~ as p — > +oo. 




As a matter of fact, the mountain-pass lemma can also be used in dimension three to assert 
the existence of a branch parametrized by the speed. 
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Theorem 12 ([18]). Assume N = 3. There exists cq > such that there exists a non-constant 
finite energy solution u c to (TWc) for any < c < cq . Moreover, v presents a vortex ring whose 
diameter is ~ HlisMI^ as c _> o, and 

. 2vr 2 | ln(c)| 2 , , , 7r 2 ln(c) 2 
E(u c ) ' - Wl , and p(u c ) as c -> 0. 

Remark 9. One may conjecture that both approaches exactly yield the same solutions. 

In [9], we revisit the minimization under constraint method, and are able to construct the 
full branch of minimizers parametrized by p in both dimensions. 

Theorem 13 ([9]). Assume N = 2 andp > 0. There exists a non-constant finite energy solution 
Up G W^(R 2 ) to equation (TWc), with c = c(up), and p(up) = p, such that Up is solution to the 
minimization problem 

E(up) = E min (p) = M{E(v),v G W(R 2 ),p(v) = p}. 

Remark 10. In particular, Theorem 13 shows that there exist travelling wave solutions of 
arbitrary small energy. This suggests that scattering in the energy space is not likely to hold. 

In dimension three, the existence result is somewhat different. 

Theorem 14 ([9]). Assume N = 3. There exists some constant po > such that 

i) For any < p < po, 

E min (p) = mi{E(v),v G W(R 3 ),p(v) = p} = v^p, 
and the infimum is not achieved in W(M. S ). 

ii) For any p > po, there exists a non- constant finite energy solution n p G W-^IR 3 ) to equation 
(TWc), with c = c{up) and p(up) = p. Moreover, E(up ) = -E'min(po) = V^Po, and for any 
P > Po, 

E(up) = E min (p) = ini{E(v),v G W(R 3 ),p(v) = p} < V2p. 

Besides the existence of minimizers, our analysis yields also properties of the curve E m [ n as 
well as of the speed c(u p ). More precisely, we have 

Theorem 15 ([9]). i) Let N = 2 or N = 3. For any p, q > 0, we have the inequality 

\E min (p) - E mhl (q)\ < V2\p - q|, 

i.e. the real-valued function p \— > E m i n (p) is Lipschitz, with Lipschitz's constant y/2. Moreover, 
it is positive and non- decreasing on R+. 

ii) Assume N = 2. Then, the function p i— > i£ m j n (p) is strictly concave. 
Hi) Assume N = 3. Then, 

(p) = V2p, V0 < p < po, 
whereas E m \ n is strictly concave on (po,+oo). 

In particular, the function E m [ n is differentiate except possibly for a countable set of values. 
Its derivative at the points of differentiability is given by the speed c(u p ), which satisfies for any 
P>0, 

< c(u p ) < y/2. 
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It remains an open problem to determine whether the curve E m \ n is differentiable or not. This 
question is related to the problem of uniqueness, up to symmetries, for the minimizer Up, which 
is completely open as well. As a matter of fact, uniqueness for any p > of the minimizer 
would lead to the differentiability of the full curve. We believe that, if at some point E m i n were 
not differentiable, then there are at least two different minimizers with different speeds. In that 
case, the function p i— > c(up) is not single valued. However, we can prove that it is a decreasing 
(possibly multivalued) function. 

In dimension two, the function -Emm has the following graph: 



E = V2p 





Our results are in full agreement with the corresponding figure given in [46]. In dimension 
three, the numerical value found in [46] for po is close to 80. Jones and Roberts have also shown 
in [46], mainly by numerical means, that in dimension three, the branch of solutions can be 
extended past the curve E = \/2p. Its representation in the E-p diagram bifurcates at some 
point, then is asymptotic to the curve E = \^2p (see the curve E up on the diagram above). 
At this stage, there is no mathematical proof of the existence of the upper branch of solutions. 
However, it is proved in [9] that the slope of the curve at the point (po,\/2po) is strictly less 
than y/2. This leaves some hope to use an implicit function theorem to construct the curve E up , 
at least near (po,\/2po)- 

Remark 11. Jones, Putterman and Roberts [46, 45] conjectured the existence of some momen- 
tum pi such that the minimizer Up has vortices for p > pi, and has no vortex otherwise. The 
numerical value found in [46] for pi is close to 75 in dimension three. There is no evidence of 
the existence of such a number pi. In particular, we do not know if pi > po in dimension three. 
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Taking advantage of the analyticity property of finite energy solutions to (TWc), we may 
derive the following additional property of the minimizing solutions we obtained in Theorems 
13 and 14. 

Theorem 16 ([9]). Let N = 2 or N = 3, p > and assume that E m [ n (p) is achieved by Up. 
Then Up is, up to a translation, axisymmetric around axis x\. More precisely, there exists a 
function u p : R x R + such that 

Up{x) = Up(xi, \x±\), Vx = (xi,xj_) E M N . 



4.2 The (KP I) transonic limit in dimension two 

In [46, 45], it is formally shown that, if u c is a two-dimensional solution to (TWc), then, after a 
suitable rescaling, the function 1 — |u c | 2 converges, as the speed c converges to to a solitary 
wave solution to the Kadomtsev-Petviashvili equation (KP I), which writes 

d t u + ud lU + dlu-d^ l {dlu) = 0. (KP I) 

Notice that the first terms correspond to the Korteweg-de Vries equation, whereas the last term 
is a transverse perturbation. As (GP), equation (KP I) is hamiltonian, with Hamiltonian given 
by 

E KP {u)= l -f {d x uf + \! ( d ^(d2u)) 2 -U n 3 , 

and the L 2 -norm of u is conserved as well. Solitary wave solutions u(x,t) = w(x\ — at,X2) 
may be obtained in dimension two minimizing the Hamiltonian, keeping the L 2 -norm fixed (see 
[22, 23]). The equation for the profile w of a solitary wave of speed a = 1 is given by 

diw-wdiw-dfw + d^\diw) = 0. (4.1) 

In contrast with (TWc), the range of speeds is the positive axis. Indeed, for any given a > 0, a 
solitary wave w a of speed a is deduced from a solution w to (4.1) by the scaling 

w a (xi,X2) = aw(^xi,ax 2 ). 

We term ground state, a solitary wave w which minimizes the action S given by 

S(v) = E KP {v) + lf 



among all the solutions to (4.1) (see [24] for more details). In dimension two, it is shown in [22] 
that w is a ground state if and only if it minimizes the Hamiltonian Ekp keeping the L 2 -norm 
fixed. We denote Skp, the action S(w) of the ground states w. 

The correspondence between (TWc) and (4.1) is given as follows. Setting e = \/2 — c 2 and 
7] c = 1 — |n c | 2 , and performing the change of variables 

AT f \ 6 f Xl 

it is shown that N e approximatively solves (4.1) as c converges to y/2. The minimizing branch 
constructed in Theorem 13 contains transonic solutions in the limit p — > 0, as shown in the 
following proposition. 
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Proposition 3. Assume N = 2. There exist some constants pi > ; Kq and K\ such that we 
have the asymptotic behaviours 

48^2 



&KP 



p 3 - K p 4 <V2p- E min (p) < Ktf 3 , VO < p < Pl (4.2) 



Let Up be as in Theorem 13. Then, there exist some constants p2 > 0, K 2 > and K3 such that 

K 2 p 2 <V2- c(n p ) < K 3 p 2 , VO < p < p 2 . (4.3) 
We next consider the map 

N p = N £p , 

where e p = e(up) = y/2 — c(up) 2 . In [8], we prove 

Theorem 17 ([8]). There exists a subsequence (p n )neN tending to 0, as n — > +00, a ground 
state w of (KP I), and a universal constant 70 > such that, for any < 7 < 70, we have 

N Pn -> w in C°' 7 (M 2 ), as n -> +00. 

Remark 12. There is an explicit solitary wave solution to (KP I) of speed 1, namely the so- 
called "lump" solution, which is written as 

w i (x 1 ,x 2 ) = 24- 



"(3 + xf+r C 2)2- 

It is usually conjectured that the "lump" solution is a ground state. It is also conjectured that 
the ground state solution is unique, up to the invariances of the problem. If this is the case, 
then the full family (N p ) p> Q converges to wg, as p — ► 0. 

Remark 13. If u c is a solution to (TWc) in dimension three, then it is also formally shown in 
[44, 46, 45], that the function w c defined by 



w c (x) = 



xi \Flx2 a/2^3\ 2 



converges, as the parameter e = y2 — c 2 converges to 0, to a solitary wave solution w to the 
three-dimensional Kadomtsev-Petviashvili equation (KP I), which writes 

d t u + udm + d'fu - d^\diu + Sfu) = 0. 

In particular, the equation for the solitary wave w is now written as 

diw - wdiw - dfw + d^ l {dlw + Sftu) = 0. 

So far, there is no rigorous proof of the existence of a branch of solutions in the transonic limit. 
This branch of solution is conjectured in [46], and represented on our graph above as the upper 
branch E up . 



5 Related problems 
5.1 Infinite energy solutions 

So far, we have restricted ourselves to finite energy solutions. Solutions with infinite energy are 
also of interest. In dimension two, a typical example is the stationary symmetric vortices for the 
Ginzburg-Landau equation (see e.g [6]). In dimension three, for small speeds, Chiron [19, 20] 
established the existence of solutions having vorticity concentrated on helices. Such solutions 
are reminiscent to some flows for the incompressible Euler equation. 
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5.2 Nonlinear Schrodinger flow past an obstacle 



As mentioned in the introduction, the Gross-Pitaevskii equation provides a model describing 
superfluidity in the Hartree approximation. When impurities or an obstacle are present in the 
superfluid, one may model these new elements adding an external potential V to the equation. 
As done by Hakim [43] in the one-dimensional case, the flow of a superfluid past an obstacle 
moving with constant positive speed c in direction xi, may be described by the Gross-Pitaevskii 
equation with a coupling with an additional potential V, namely 

id t V(x,t) = A^(x,t)+^(x,t)(l-\^(x,t)\ 2 -V(x 1 -ct,x_ L )\ (x = (x u x ± ),t) EK^xl. (5.1) 

As before, one adds the condition |^(a;)| — ► 1, as \x\ — ► +oo, to describe a superfluid which is 
at rest at infinity. In the frame of the moving obstacle, equation (5.1) may be recast as 

id t $ = icdi<$> + A$ + $(l - |$| 2 - v\ (5.2) 

where we denote 3>(cc,i) = *$>(x\ + ct,x±,t). Stationary solutions v to (5.2) satisfy the elliptic 
equation 

icd!V + Av + v(l -\v\ 2 -V^J =0. (5.3) 

Equation (5.2) formally remains hamiltonian. For a superfluid at rest at infinity, the conserved 
Hamiltonian may be written as 

F?m = \l \v<s>\ 2 + \[ (i-m 2 ) 2 -U v{i-\^\ 2 )- c -f ^,^-1). (5.4) 

2 Jrn 4 J r n 2 J r n 2 J r n 

The previous variational formulation leads us to restrict ourselves to solutions v for which the 
following modified Ginzburg-Landau energy is finite, that is 

E V {v)^\f \Vv\ 2 + \[ (1 - \v\ 2 ) 2 - \ [ y(l-H 2 )<+oo. 
2 Jrn 4 J r n 2 J r n 

In [43], a formal and numerical analysis of stationary solutions to (5.2) was initiated in the 
one-dimensional case, whereas in [54], still in dimension one, Maris established the existence 
of solutions to (5.3) minimizing the Hamiltonian FY , assuming the potential V is a bounded 
measure with small total variation compared to the speed c. Assuming additionally that V is 
compactly supported, Maris [54] also exhibited excited states. 

In higher dimensions, the flow of a superfluid past an obstacle has been widely investigated 
both formally and numerically (see, for instance, [27, 59, 60, 62, 1, 2, 16]). In [64], Tarquini 
proved that, under suitable assumptions on the localized potential V, there are no finite energy 
supersonic solutions to (5.3), that is with speed c > y/2. He also computed an explicit bound 
for the solutions to (5.3), and described their regularity and convergence at infinity. 

To our knowledge, the existence problem of finite energy solutions to (5.3) has not been 
addressed yet for a general class of potentials, although some special cases have been considered 
in [4, 3]. Therefore, we would like to provide in this paper an existence result in dimension two. 

Theorem 18. Let < c < y/2 and V S L 2 (K 2 ). There exists a constant K(c) > depending 
only on c, such that, if 

[ \V\ 2 <K(c) 2 , (5.5) 
then equation (5.3) has at least one finite energy solution u c . 
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The solutions u c provided by Theorem 18 are constructed minimizing the Hamiltonian FY 
given by (5.4) locally near the constant solutions to (TWc). As a matter of fact, in the case 
the potential is 0, the solutions we construct correspond to the constant solutions to (TWc). So 
far, we have not tried to construct excited states (as done in the one-dimensional case by Maris 
[54]). For a zero potential, these solutions would correspond to the non-constant two-dimensional 
solutions to (TWc) we have discussed before. 

Remark 14. Though we did not address this question here, we believe that the solutions u c 
provided by Theorem 18 are local minimizers of FY on suitable variational spaces (see [9], where 
a similar question is addressed for the equation without potential, and [64] for the choice of the 
suitable variational space, in case V has compact support). In particular, the argument used to 
prove Theorem 18 does not establish that either FY(u c ) or p(u c ) are finite. 

As mentioned above, we construct the solutions u c by locally minimizing the Hamiltonian 
FY given by (5.4). When the potential V belongs to L 2 (M. 2 ), FY is well-defined and of class 
C 1 on the space {1} + if 1 (IR 2 ). However, the construction of minimizers for FY presents some 
difficulty in view of the lack of compactness. To circumvent this difficulty, we introduce, as in 
[10, 9], minimizing problems on expanding tori, for which the existence of minimizers presents 
no major difficulty, and then pass to the limit when the size of the torus tends to infinity. More 
precisely, we introduce the flat torus defined by 

T n ~ 0„ = [-• 7rn, 7rn] 2 , 

for any n £ N* (with opposite faces identified), and the space 

X n = H 1 (T n ,C)-Hl CT (n n ,C) 

of 2-7rn-periodic .ff ^functions. We define the energy EY and the momentum p n on X n by 

E V n(v) = \f Wv\ 2 + \[ (l-\v\ 2 f-U V(l-\v\% (5.6) 

and 

Pniy) = - {idtv,v- 1), 
z Jn n 

so that the Hamiltonian FY n on X n is given by 

FY n (v) = EY(v)-c Pn (v). (5.7) 

Notice that the function V in (5.6) and (5.7) denotes the 2irn -periodic restriction of V to Q n . 
For given A > 0, we introduce the open set £ n (A) defined by 

1 ' — |2 , 1 / /i i..|2\2 



£ n {A) = |d£ X n , s.t. E n {v) = 2j a \ Vv \ +lj n ( X " \ y n < A /' 
which is a neighbourhood of the set of constant functions, and consider the minimization problem 



J^ B (A)= inf [FY n (v) . (5.8) 

Our purpose is to prove, that for a suitable choice of A = A(c) with respect to the speed c, 
independent of n, the previous minimization is achieved. The main ingredient in the proof of 
Theorem 18, which allows in particular to prove that FY is bounded from below on any set of 
functions with sufficiently small energy, is the following estimate. 
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Lemma 17. Let 5 > and 8\ > such that 

5 c = -j=<5 <5 1 <1, 

and consider a function v G £ n (A). Then, there exists some integer uq, and some positive 
constant Kq, depending only on do and 5\, such that 

c\ Pn (v)\ < 5 fE n {v) + — g 1 jr E n (v)l + 2 ^ C £» 2 , (5.9) 

d 0r(l - 5f )(<5i - d ) (1 - 5f){8 1 - d ) 2 d 

for any n > no, and any A < Ao- 

Proof. The proof of Lemma 17 is reminiscent of the proof of Lemma 4. When E n {v) is small 
enough, the modulus of v is close to 1. Assuming that we may construct a lifting v = gexpiip 
of v, we compute 



c\Pn{v)\ ~ | 



/ {l-Q 2 )d W <^( ^«d c / e(v)=5 c E n (v). (5.10) 
J!J„ V2 Jn n 8 Jn„ 



Lemma 17 reduces to estimate any error terms in (5.10) to obtain (5.9). 

Therefore, we consider some function v G £ n (A), assuming, up to a density argument, that 
v is smooth on O n . In order to estimate how the modulus of v is close to 1, and to construct a 
lifting of v, we study the level sets 

ft(0 = {x G fi n , s.t. =£}, 

for any given £ G R + . By Sard's lemma, the sets fi(£) are one-dimensional, smooth submanifolds, 
for almost any £ G R+, so that their Hausdorff's length 

L(o=n 1 (m), 

is finite for almost any £ G R+. Moreover, the coarea formula gives 
Si / r \ r 

2, 



|Vi/|]df< / |Vu(x)| 2 dx < 2£ , n (u) 

<5 Wn(£) / i<5 <|x|<<5i 



and 

rS 



/ [ \Vv(x)\dx< V^^nW^C/^l)^, 

J5 J(5 <|x|<5i 

where we denote [/(£) = {x G ri n , s.t. |f(x)| < £}, for any £ G R+. Since 

Wl)l < JY^ E n« 

we are led to 

f 5l L(0dt<^E n (v). 

J So l ~ d l 

The mean value inequality then provides the existence of do < £o < di such that the length of 
the one-dimensional submanifold fi(£o) verifies 

/ x , \ 2 V / 2A 

(1 - df)(di - d J (1 - df)(di - d ) 
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and 



/ 



2 , , 2A 
Vv < -Ejv) < 



8 - 5 1 



So - 6 1 ' 



(5.12) 



We first deduce from (5.11) that, up to some new unfoiding of the torus T n , f2(£o) is, for 
n > w ( 1 Jtfr^ i _ So <) , a collection of smooth, connected, closed curves (7i)ie/ of Q n . Moreover, 



since 



i^(6)i < < 



4 „ , . 4A 
E n (v) < 



(1 " *l) 2 



(i -sir 



we may assume that the complementary of U (£o) contains the boundary of fi n for any n > 

It now follows from (5.12) that 

u(y) If ,_ , . 2A 



sup 



|u(ai)| |v(y)| 



< 



iVul < 



& jQ(^o) ^o(<^l - S ) 



so that, assuming that A < <5o(<5i — #o), the topological degree of A on each curve 7$ is equal 
to 0. This provides the existence of a smooth lifting v = gexpip 01 u on the complementary of 
U(£o), so that, using the periodicity of v, we can compute 



1 



Pn(v) = - I (idiv,v) = - I (idiv,v)-- 



1 



g 2 dnp. 



n n \u(£ ) 



(5.13) 



The first term in the right-hand side of (5.13) can be estimated by 



Ufa) 



(id±v, v) 



< 



C/(€o) 



so that, using the isoperimetric inequality 

47r|[/(£o)| < L(e ) 2 , 

and (5.12), we are led to 



Ufa) 



(idiv,v) 



< 



Si 



V5F(1- Wi-*o) 



(5.14) 



On the other hand, the second term in the right-hand side of (5.13) is written as 



If Q 2 dw=\ I (l-Q 2 )dnp- -J2 I (Y>-<Pi)vi, ( 5 - 15 ) 

1 Jn n \u(£ ) 1 Jo n \u(to) 1 ieI Ji, 



where v\ denotes the first component of the outward normal to 7$, and 

1 



Using the result of Lemma 4, we obtain 



n„\t/(£o) 



(1 - Q 2 )d W 



< V2 



ip. 



e(v) , . 

n„\c/(So) 8 £0 



whereas 



£ / (v-^i <£ / ( I i v ^i) 



< 



Hto 

£0 Jn(to) 



\Vv\ 
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so that using (5.11), (5.12) and (5.15), we are led to 



finW(£o) 



< 



V25 { 



■E n (v) + 



2^2 



*o(l-<5?)(*i-*o) a 



EJvY. 



(5.16) 



Inequality (5.9) finally follows from (5.13), (5.14) and (5.16), provided that 



A < Aq = Sq(8i — 5q), and n > no = max 



As a consequence of Lemma 17, we have 



2V25 2V<5 ($i-(So) 



vr(l-Jf)' n(l-8() 



(5.17) 
□ 



Lemma 18. Let < c < \/2 and n > no, where no is i/ie integer of Lemma 17. There exists 
positive constants K{c), A(c) and n(c) < A(c) such that, if \\V\lL2m2\ < K(c), then 



<„W > ^„(«(c)), 
/or any function v £ £ n (A(c)). Moreover, 

-A(c) < ^ c y n(A(c)) < 0. 



(5.18) 



(5.19) 



Proof. The proof reduces to compute an upper and a lower bound for the function F^ n given 
by 

(5.20) 



,)=E n {v)- 1 - t V(l- 

J T n 



on the set £ n (A), which amounts to estimate the second and the third term of the right-hand 
side of (5.20). We use Lemma 17 to provide an upper bound for the scalar momentum p n (v). 
Therefore, we consider numbers 6q, 5i, 82 and // such that 



c 5 C 
5 C = —= < 82 < 5q < 8\ < 1, and < fi < 1, 
V2 o 2 



(5.21) 



and assume that n > no and A < Aq, where no and Aq are provided by Lemma 17. By Lemma 
17, we have 

\cp n (v)\ < Y-E n (v), 
82 

provided that E n (v) < Ai(c), where Ai(c) denotes some positive constant depending only on c. 
On the other hand, we compute 



V(l - \v\ 2 ) 



< 



4( 



r^)I/ 2 + (1 -^ 



so that 



^|)^ (U) -4(1^) 



We now fix 



/ V 2 <F^ n (v)< (l + ll+ 5 f)E n {v) + — — [ V 2 . (5.22) 
Jm? v 82' 4(l-/i)7 R2 



A(c) = min{A , Ai(c)} > 0, 
and choose some positive number K(c) and k(c) such that 



(1 + 2h)k(c) < A(c), and K(c) 2 < 4(1 - fi) U - ^)«(c), 



(5.23) 
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so that, assuming that 1 1 ^ 1 1 jl 2 (m 2 ) — K{ c )i an d n — n 0i (5.22) provides 

F c » > (/i - |)k(c) - i ^yA'(c) 2 > = Fl n (l) > ^ n (K(c)), 

for any function i> G AT n such that k(c) < E n {v) < A(c). Hence, we are led to (5.18), that is 

^k(c))=^(A(c))<0. (5.24) 
On the other hand, (5.22) also yields, using (5.23), 

\Fl n {v)\ < (l + + |) «(c) + 4(1 1 _ /i) ^(c) 2 < (1 + 2p)s(c) < A(c), 

for any function u G £„(k(c)), so that (5.19) follows from (5.24). □ 

We next have 

Lemma 19. Let < c < \pl and n > hq, where uq is the integer of Lemma 17, and let K(c) and 
A(c) be as in Lemma 18. Assume moreover that condition (5.5) holds for K(c) and V . Then, 
the minimization problem (5.8) for A = A(c) has a solution u™ in £ n (A(c)). In particular, u" is 
solution to (5.3) on ft n , and verifies 

E n (u n c ) < k{c) < A(c), and - A(c) < F^ n (u n c ) = ^ c y n (A(c)) < 0. (5.25) 

Proof. Lemma 19 follows from standard variational arguments. For given n G N and A > 0, 
F^ n is bounded from below on £ n (A). Indeed, there exists some constant K(n, A) such that 

/ |V?;| 2 + I \v\ 4 < K(n, A), 

for any function v G £ n (A), so that 

F c y (v) > -VX(\\V\\ L 2 { j n) + c\\v-l\\ L 2 {Jn ^ > —K(n, A). 

Next, consider a minimizing sequence (w m ) m ^ for i* 1 ^ on A(c), and assume n > no- In view 
of (5.18), we have 

?nW > •^(^(c)), 

so that we may assume, since < k(c) < A(c), and hence £ n (n(c)) C £ n (A(c)), that u; m G 
£n(^(c)) for any m G N, or equivalently, 

KK) < k(c) < A(c). 

Using Rellich's compactness theorem, there exists, up to some subsequence, a function u™ G X n 
such that 

Vw m — 1 Vzi™ in L 2 (Q n ), and w m — > u™ in L 4 (Q n ), asm-* +oo, 

so that 

Fl n (u n c ) < liminf fe( Wm )) = ^ n (A). 

The same argument also provides 

£ n «) < liminf (E n (w m )) < n(c) < A(c). 

when ||U||^2(jg2) < K(c). Hence, the minimization problem (5.8) for A = A(c) is attained by it™ 
in the open set £ n (A(c)). Since F^ n is of class C 1 on X n , u™ is a critical point of that is a 
solution to (5.3) on Q n . □ 
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Proof of Theorem 18. Passing to the limit n — > +00 requires to prove some compactness for the 
sequence (ii") n > no . We first provide a uniform estimate for (u^) n > no using (5.25), and standard 
elliptic estimates. By (5.25), 



[ (|V<| 2 + K| 4 ) <A(c) 



'B(x ,2) 

where we denote B(xo, 2) = {x G Q n , \x — xq\ < 2} for any xq £ Q n . Using Holder's inequality, 
we are led to 

/ (|<(i - ki 2 )|^ + + \icd lU ?fA < a(c), 

JB(xo,2) V ' 

so that, by equation (5.3), and standard elliptic estimates, the sequence (u™) n >n is bounded in 
W 's(B(xo, !))• Therefore, by Sobolev's embedding theorem, there exists some positive constant 
A(c) such that 

IKIIloc^) <A(c), Vn>l. (5.26) 

Moreover, (5.25) also gives 

||V<|| i2(nn) < A(c), Vn > 1. (5.27) 

Since any compact set K of M? is included in O n for n sufficiently large, we then construct, 
using (5.26), (5.27), Rellich's compactness theorem, and a diagonal argument, a function u c £ 
H\ oc (M. 2 ), such that, up to some subsequence, 

Vu™ Vu c in L 2 (K), and u n c -> u c in L 9 (liQ, asm +00, (5.28) 



for any compact set A" of M 2 , and any 1 < q < +00. Since is solution to (5.3) on fi n , 
convergences (5.28) first yield that u c is solution to (5.3) on M 2 . Moreover, they also provide 



2 

so that by monotone convergence 



#K) < A(c). 



Hence, 

^ y K) < £(u c ) + \\V\\ L 2 m E(u c )^ < +00, 
so that u c is a finite energy solution to (5.3) onf 2 . □ 

Remark 15. Notice that it follows from (5.17), (5.21) and (5.23) that the constant K(c) in 
Theorem 18 tends to 0, as c — > so that the considered potential V must be smaller and 
smaller as c — > y/2. 



5.3 Nucleation by impurities 

An interesting generalization of (GP) (see [32]) concerns the modelling of vortex nucleation by 
an impurity, e.g. an electron. In the Hartree approximation, the equations governing the one- 
particle wave function of the condensate ip, and the wave function of the impurity <fi, are a pair 
of coupled equations 

h 2 

ihd t V = -7H7 A ^ + (^o|$| 2 + ^ol^l 2 - E)V, 

2 f (5.29) 

ihd t <& = -— a$ + (f/ |*| 2 - E e )$, 
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where M and E are the mass and single-particle energy for the bosons, and (i and E e are the 
mass and energy of the impurity, while Uq (resp. Vo) denotes the mass of the Dirac interaction 
potentials between bosons and impurity (resp. bosons). 
In dimensionless variables, (5.29) is written as (see [32]) 



id t v = -a* + 1(1|$| 2 + 1^| 2 - 



i5d t <5> = -A$ + 4j (V|^| 2 -e 2 k 2 ^j<5>. 



(5.30) 



Here, 8 = is the ratio of the mass of the impurity over the boson mass (5 <C 1), q 2 = 

k 10 = 2n 2 M%E$h ,i 1S a dimensionless measure of the single-particle impurity energy, and e 10 = 

emHj" 1 1S a dimensionless constant, which in applications is about 0.2. 

Assuming that we are in a frame in which the condensate is at rest at infinity yields the formal 
boundary conditions 

|V>(aO| — ► 1, and <f>(x) — ► 0, as \x\ — > +oo. 
Similarly to the Gross-Pitaevskii equation, (5.30) has a formally conserved Hamiltonian, 

Iff E 2 e 2 / \ 2 e 2 k 2 \ f 

E(V,$) = - f e 4 |V^| 2 + ^r|V4>| 2 H 1 - l^l 2 + |^| 2 |*| 2 H$| 2 =/ e(*,$). 

2 Jrn V q l 2 V J q 2 J Jrn 

It moreover conserves 

|2 



m($) = / |$| 



The finite energy travelling waves of (5.30) have been formally investigated in [32]. No rigorous 
existence results seem to be known so far in dimensions two and three, but Bouchel [15], ex- 
tending the analysis of [34, 35] for (GP), has proved decay estimates for finite energy travelling 
waves of (5.30), together with the non-existence of supersonic travelling waves. On the other 
hand, Maris [55] has proved that one-dimensional travelling waves to (5.30) exist if and only if 
their velocity is less than the sound velocity at infinity, and that in this case, the set of travelling 
waves contains global subcontinua in appropriate Sobolev spaces. 



6 Conclusion 

Although travelling waves have been widely discussed in the physical literature, the field is 
still widely open for mathematical investigation. In the course of our discussion, we already 
mentioned a number of open problems and research directions: differentiability of the curve 
E m - m , existence of solutions to (TWc) for all values of speeds < c < y/2, uniqueness, up 
to the invariances, of constrained minimizers as well as of ground states of (KP I), existence 
of the upper branch E up in dimension three, and its possible transonic limit, non-existence of 
sonic travelling waves in dimension three, similar questions for the related problems above... To 
complete this survey, we wish to emphasize the question of stability: even the orbital stability 
of constrained minimizers presents serious yet unsolved difficulties in dimensions two and three. 
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